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ABSTRACT 

Gravitational lensing by clusters of galaxies has been detected on scales 
ranging from ~ 1CT 1 Mpc to ~ 10 Mpc, namely, arcs/arclets, weak lensing and 
quasar-cluster associations. This allows us to derive an overall radius matter 
distribution of clusters of galaxies. While the dynamical analysis of the X-ray 
observations has yielded a great number of data for the virial cluster masses, it 
becomes possible to statistically compare the cluster mass determinations by 
these two independent methods. In this letter we show that as compared with 
gravitational lensing, the dynamical analysis under the assumption of isothermal 
and hydrostatic equilibrium has systematically underestimated the cluster 
masses inside the Abell radius by a factor of ~ 2 with scatter between 0.7 and 
5. Because the same correction factor should be applicable to the gas baryon 
fraction of clusters of galaxies obtained from the X-ray data, it is probably too 
premature to claim a baryon crisis in today's cosmology. 

Subject headings: cosmology: theory — galaxies: clusters: general - 
gravitational lensing 
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1. Introduction 

A combination of the primordial nucleon abundances predicted by the standard Big 
Bang Nucleosynthesis (BBN) and those inferred from astronomical observations has set a 
tight constraint on the baryonic matter component of the universe (Walker et al. 1991): 
0.04 < f^/iJo < 0.06. This indicates that the baryon fraction, f b = £l b /£l with f2 being the 
average mass density of the universe in units of critical density, is smaller than ~ 0.06 
in the prevailing cosmological model of Q = 1, and a significant fraction of the mass in the 
universe should be invisible (non-baryonic matter). However, such a standard scenario has 
been challenged in recent years by the X-ray observations which detect a considerably large 
amount of the hot X-ray emitting gas in clusters/groups of galaxies. The resulting gas 
baryon fraction is a few time greater than the prediction of BBN, provided that the gas is 
in the state of hydrostatic equilibrium with the gravitational potentials of clusters/groups 
of galaxies. The baryon crisis thus arises if the matter in clusters/groups of galaxies is 
representative of the universe. In particular, this discrepancy probably implies that at 
least one of the basic hypotheses in our current theories of cosmological study needs to be 
modified or even abandoned (White et al. 1993). 

Yet, the above claim should be taken very cautiously without carefully examining the 
reliability of the X-ray cluster mass determinations. Indeed, the existence of substructures 
and the recent detection of the complex two-dimensional temperature patterns in clusters 
of galaxies (e.g. Henry & Briel 1995; Markevitch 1996; Henriksen & Markevitch 1996; 
Henriksen & White 1996) strongly suggest that clusters of galaxies may not be the 
well-virialized dynamical systems as were believed before and the uncertainty in cluster 
mass determinations assuming hydrostatic equilibrium for the X-ray gas may be quite large 
(Balland & Blanchard 1996). Therefore, it is desirable that another independent cluster 
mass estimate is made to test the accuracy of the X-ray cluster mass determinations and 
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furthermore, to re-examine whether there is a baryon overdensity in clusters of galaxies. 

It has been realized that gravitational lensing associated with clusters of galaxies 
can fulfill the task, which gives rise to cluster masses regardless of the cluster matter 
components and states. In several clusters of galaxies where both X-ray data and image 
distortions of background galaxies are available, comparisons of the virial cluster masses 
derived from X-ray observations and the gravitational cluster masses inferred from the 
distorted images of distant galaxies have been made (e.g. Wu 1994; Fahlman et al. 1994; 
Miralda-Escude & Babul 1995; Squires et al. 1995). Today, gravitational lensing by clusters 
of galaxies has been detected on scales ranging from the inner core to the outer radius 
of ten arcminutes, including giant arcs/arclets, weakly distorted images of background 
galaxies and quasar-cluster associations. These lensing observations alone may allow us to 
derive an overall radius matter distribution of clusters. Therefore, it would become possible 
to statistically compare the cluster matter distributions given by dynamical analysis of 
the X-ray observations with those by gravitational lensing. This procedure is essentially 
different from the previous work which focused on individual clusters with both X-ray and 
gravitational lensing observations. We now select the two sets of data separately from 
literature. This letter presents the result of the comparisons and discusses its significance for 
cosmological study. Throughout this letter we adopt a matter-dominated flat cosmological 
model of Q = 1 and a Hubble constant of H = 50 h 50 km s _1 Mpc. 

2. X-ray cluster mass determination and gas baryon fraction 

Under the assumption of the standard isothermal /5-model for the X-ray surface 
brightness of cluster of galaxies, the total mass in gas within radius r of cluster center 
is (Cowie, Henriksen, & Mushotzky 1987) M gas (r) = Annor^fimp 

while the equations of hydrostatic and dynamical equilibrium give the total virial mass 
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M vir (r) = 3/3(kT/ [im p G)r 3 / (r 2 +r 2 ), where n and r c denote, respectively, the central 
number density and core radius of the gas profile, T is the gas temperature, k, the 
Boltzmann's constant and [im p , the mean particle mass. The ratio of M gas to M V i r provides 
a conservative estimate of the cluster baryon fraction /& since the galaxy contribution is not 
included. Fig.l shows the measured gas baryon fractions of clusters of galaxies in literatures 
without any corrections, in which we have only utilized the virial masses obtained in the 
case of isothermality. Note that for most of the clusters the gas baryon fractions at radii 
of larger than ~ 1 Mpc are computed from the spatially-unresolved measurements of the 
gas temperature. This leads to an underestimate of gas baryon fraction if temperature 
decreases with radius as it is naturally expected (e.g. Henriksen & Mamon 1994). The 
relatively low gas baryon fraction at the largest radius r ~ 7.1 - 10 Mpc for A2142 in 
Fig.l probably arises from such an oversimple assumption (Henriksen & White 1996). It 
appears that based on the current data, we have not detected any apparent variations of 
the gas baryon faction of clusters of galaxies with radius. The mean gas baryon fractions 
are f b « 14% and f b ~ 18% with and without those data at r = 0.5 Mpc given by Edge & 
Stewart (1991). It is thus concluded that the mean baryon fraction in clusters of galaxies is 
about 2-4 times larger than the prediction of BBN, if the hot X-ray gas is in hydrostatic, 
isothermal equilibrium with the binding cluster gravitational potentials. 

EDITOR: PLACE FIGURE HERE. 



3. Cluster masses from gravitational lensing 



Arcs/arclets are the strongly/moderately distorted images of distant galaxies by 
foreground clusters of galaxies. The projected cluster mass within the position (r arc ) of 
arc/arclet can be easily obtained if one assumes a spherical matter distribution for the 
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lensing cluster and approximates the alignment parameter of the background galaxy to zero: 
miens(r arc ) = 7trl rc j: crit , where T, crit = (c 2 /AtiG){D s / D d D ds ) is the critical mass density with 
Dd, D s and Dd s being the angular diameter distances to the cluster, to the galaxy and from 
the cluster to the galaxy, respectively. For the complex arc/arclet configurations, cluster 
mass can be estimated by constructing an asymmetrical lens model [see Fort & Mellier 
(1994) for a recent review]. We illustrate in Fig. 2 the cluster masses given by modeling of 
arcs/arclets. One major uncertainty comes from the unknown redshifts for some arclike 
images, for which we have assumed z arc = 0.8. For a typical arc-cluster at redshift of ~ 0.3, 
this leads to an overestimate of cluster mass by a factor of 1.4 if the background galaxy is 
actually located at z arc = 2. 

Another powerful tool of probing the matter distribution of cluster of galaxies is to 
study the weak gravitationally induced distortions in the images of faint galaxies behind 
cluster of galaxies. By analyzing the shear field (tt) around the cluster, the statistics 
C( r ) — Ir max (lT)(^ — r2 / r marr) _1 ^ mr measures the mean surface mass density in units 
of E cr j t interior to r minus that in the annulus from r to r max (Fahlman et al. 1994). 
Therefore, a lower bound on the projected cluster mass within the radius r can be found 
through mi ens {r) = 7rr 2 £(r)£ cr jf. Cluster mass reconstructions have been made for several 
clusters of galaxies in which the statistically significant shear patterns are detected. The 
resulting cluster masses are plotted in Fig.2. Again, there has been so far no information 
available about the redshifts of background galaxies and a mean value of (z) — 1 - 3 has 
been often assumed in the computations. This brings about an uncertainty of cluster mass 
by a factor of ~ 1.3 for a lensing cluster at redshift of ~ 0.3. 

Gravitational magnification can also enhance the number density of background sources 
around a foreground cluster of galaxies, which has been recently confirmed by discovering 
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the so-called quasar-cluster associations on scale of up to ~ 10 arcminutes (Wu & Fang 
1996 and references therein). One can figure out the mean cluster mass which is required 
to produce the reported quasar overdensity in terms of gravitational lensing. It turns out 
that clusters of galaxies should contain considerably large gravitational masses extending 
to a radius of ~ 10 Mpc in order to account for the quasar enhancements. We compute the 
projected cluster masses over the association areas simply by mi ens (r) = -nr 2s L and show the 
results in Fig. 2, where the mean cluster surface mass density £ have been given by Wu & 
Fang (1996) for the four measurements of quasar-cluster associations. 

It is remarkable that the projected cluster masses revealed statistically by three 
different lensing methods over two decades in radius from ~ 10 _1 Mpc to ~ 10 Mpc 
can be well fitted by a power-law: mi ens (r) = l0 15 - 39±0 - 17 (r/Mpc) L51 M , where (also 
hereafter) the error bar represents the scatter of the best-fit average value rather than 
the real uncertainty in the measurement which is difficult to estimate. Because the weak 
lensing method usually provides a low limit to the cluster mass, most of its results are 
smaller than the mean value. The fitting without weaking lensing data becomes slightly 
steeper: mi ens (r) = io 15 - 56±0 - 11 (r/Mpc) 1 - 63 M . While one might argue the reliability of the 
cluster masses up to the radius of r « 10 Mpc derived from the quasar-cluster associations, 
we give the fit by removing the four results of quasar-cluster associations from Fig. 2: 
m lens (r) = 10 15 - 22±ai5 (r/Mpc) L32 M . 



4. Comparisons 

We also display in Fig. 2 the X-ray cluster masses M vir (r) derived from the isothermal 
/9-model under the assumption of hydrostatic equilibrium, including the results for five 
compact groups of galaxies. Recall that M vir (r) are the dynamical masses used in Fig.l for 
computations of the gas baryon fractions in clusters of galaxies. The discrepancy of M vir 
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and mi ens (r) is clearly seen at small radius, while the two sets of data seem to merge beyond 
r ~ 1 Mpc. Therefore, an intuitive speculation for such a variation is the projection effect. 
We have then tested the conventional r~ 2 profile and the so-called universal density profile 
found by the standard CDM simulations (Navarro, Frenk & While 1995). Neither of these 
profiles can fit both the three-dimensional masses and the two-dimensional projected ones. 
In fact, the deprojection of rrii ens (r) recovers the corresponding three-dimensional masses 
Mi ens (r), if one assumes a spherical matter distribution as we have already adopted in the 
above sections. We list the resulting Mi ens (r) in Table 1, together with a least-square fit of 
a power-law to the X-ray cluster masses M vir {r) and the ratios of Mi ens {r) / M vir at different 
cluster radii. 



EDITOR: PLACE FIGURE | HERE. 



EDITOR: PLACE TABLE g HERE. 



There is a significant disprepancy between the virial and lensing cluster masses inside 
the core radius of the X-ray gas profile, r c ~ 0.25 Mpc, where arcs/arclets are observed. It 
seems that lensing method using arclike images can give rise to the cluster masses of ~ 5 
times larger than the virial masses. This is essentially comparable to the previous similar 
studies for individual clusters (Wu 1994, Miralda-Escude & Babul 1995). Around the radius 
of r ~ 1 Mpc, at which the lensing data are dominated by the weak lensing observations, 
the mean ratio of Mi ens /M V i r is about 2-3. However, considering the large scatters, one 
cannot exclude the possibility that the cluster masses derived from the two methods are 
consistent. Yet, most of the data based on the weaking lensing method correspond to the 
low limits to cluster masses. The similar conclusion can be applicable to the Abell radius of 
3 Mpc, where one may expect a mean factor of ~ 2 in Mi ens /M V i r . The mass discrepancy 
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might vanish at the outer radii of clusters (r > 5 Mpc) if we take out the results given 
by the quasar-cluster associations which dominate the cluster mass determinations with 
gravitational lensing at large cluster radius. 

5. Discussion and conclusions 

It appears that as compared with the gravitational lensing method, the dynamical 
analysis based on the isothermal, hydrostatic equilibrium has systematically underestimated 
the cluster masses within the Abell radius (3 Mpc) by a factor of ~ 2 with scatter ranging 
from 0.7 to 5. The mass discrepancy is rather remarkable inside the cluster core of r c ~ 0.25 
Mpc but diminishes along the outgoing radius. As an immediate result of this discrepancy, 
the baryon fractions in clusters of galaxies provided by the X-ray cluster masses should 
be correspondingly reduced by the same factor, which thus opens a possibility to solve or 
partially remove the recently claimed baryon crisis in clusters of galaxies. Meanwhile, our 
finding indicates that clusters of galaxies may not be regarded as the well-relaxed virialized 
systems. 

The above conclusions are strongly supported by the spatially resolved spectra for 
some clusters of galaxies obtained with ASCA, GINGA and ROSAT, which show the 
complex temperature patterns over the cluster faces (see Henry & Briel 1995; Markevitch 
1996; Henriksen & Markevitch 1996; Henriksen & White 1996). These significant 
temperature variations cannot be described by a simple analytic profile like a /3-model, 
and non-isothermality in the hot X-ray emitting gas of clusters of galaxies is apparently 
required. So, cluster mass determinations using the isothermal hypothesis for the X-ray gas 
may lead to large errors. 

Gravitational lensing is a robust estimate of gravitational mass in a celestial body 
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or system. We have found the consistency between the cluster masses derived from 
three different lensing phenomena over scale 0.1 Mpc < r < 10 Mpc. The total masses 
inferred from lensing can be well represented by a single power-law of ~ r i-5±o.2^ i nc [i ca ti ve 
of a density profile of ~ r -i-5±o.2 rp^jg - g s ^ ee p er than the singular isothermal matter 
distribution ~ r~ 2 and lies between the ~ r~ l and ~ r~ 3 universal profile predicted by the 
standard CDM simulations (Navarro et al. 1995). 

Nonetheless, the mass discrepancy between dynamical analysis and gravitational 
lensing at the central regions of clusters of galaxies could arise from the cooling flow and/or 
the contribution of the nonthermal pressure such as magnetic field (Loeb & Mao 1994; 
Ensslin et al. 1996). Alternatively, gravitational lensing may overestimate cluster masses if 
the lensing cluster is prolate with the major axis along the line-of-sight (Miralda-Escude & 
Babul 1995). Furthermore, it is somewhat hard to understand the cluster mass extension 
to a radius of ~ 10 Mpc where the quasar-cluster associations are detected. 

After all, neither of the present X-ray cluster masses and gravitational lensing data 
forms a complete sample and therefore, it is impossible to evaluate the statistical significance 
of the result. What we would like to emphasize in this letter is that the comparisons of 
the updated cluster mass determinations by these two independent methods have revealed 
the possible existence of large mass uncertainties under the scenario of the hydrostatic 
equilibrium. So, it is probably too premature to claim a baryon crisis in today's cosmology. 

We thank an anonymous referee for helpful suggestions. WXP was supported by the 
National Science Foundation of China and a World-Laboratory fellowship. 
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Table 1: Comparisons of virial and lensing cluster masses. 



method 


power-law fit 


Mi ens /M vir 


r = 1 Mpc r = 3 Mpc r = 5 Mpc 



virial M vir {r) = io 14 - 53±0 - 23 ( r /Mpc) L91 M ( 



lensing (1) M lens (r) = 10 14 - 96±0 - 06 (r/Mpc) 1 - 51 M 2.69l£i$ 1.73lJ;£ 1.411J;; 
lensing (2) M ZeriS (r) = 10 14 - 90±0 - 03 (r/Mpc) 1 - 32 M 2.341^1 1.231^° 0.9118: 
lensing (3) M ieriS (r) = 10 15 - 04±a07 (r/Mpc) L63 M 3.24±f;|l 2.38l?$ 2.061?;! 



(1) Arclike images, weak lensing and quasar-cluster associations; (2)Arclike images and weak 
lensing; (3)Arclike images and quasar-cluster associations 
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Fig. 1. — Baryon fractions in clusters of galaxies derived from the X-ray observations 
under the assumption of isothermal and hydrostatic equilibrium. The dashed lines show 
the predictions of the standard nucleosynthesis for the cosmological models of Q = 1 and 
VI = 0.3. (see Fig. 2 for references) 

Fig. 2. — Comparisons of cluster masses inside radius R derived from dynamical analysis 
and gravitational lensing. Open squares: the projected cluster masses from arcs/arclets 
(Wu 1994; Kneib k Soucail 1995); Fancy squares: the projected cluster masses from weak 
lensing technique (Fahlman et al. 1994; Tyson k Fischer 1995; Smail k Dickinson 1995; 
Luppino k Kaiser 1996; Seitz et al. 1996; Squires et al. 1996a,b); Diamonds: the projected 
cluster masses from quasar-cluster associations (Wu k Fang 1996); Octagons: the cluster 
masses from X-ray observations (Hughes et al. 1989; Edge k Stewart 1991; Briel, Henry, 
k Bdhringer 1992; Miyaji et al. 1993; White et al. 1993; Briel k Henry 1994; White et al. 
1994; Elbaz, Stewart, k Bdhringer 1995; Dell'Antoio, Geller, k Fabricant 1995; White k 
Fabian 1995; Schindler 1995; Schindler k Wambsganss 1995; Ikebe et al. 1996; Bohringer 
et al. 1996; Schindler et al. 1996; Squires et al. 1996a; Henriksen k White 1996). The 
masses of compact groups of galaxies from X-ray observations are also illustrated (crosses) 
as a comparison (Henriksen k Mamon 1994; Pildis, Bregman, k Evrard 1995). The solid 
lines show the best fitting power-laws of the dynamical and lensing data, and the dashed line 
is the three-dimensional cluster masses obtained by the deprojection of the two-dimensional 
lensing masses. 
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